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CN ■ Abstract: We study how the universe reheats following an era of chaotic Dirac-Born- 

Infeld inflation, and compare the rate of particle production with that in models based on 
canonical kinetic terms. Particle production occurs through non-perturbative resonances 



in 

^D I whose structure is modified by the nonlinearities of the Dirac-Born-Infeld action. We 

^-^ ' investigate these modifications and show that the reheating process may be efficient. We 



estimate the initial temperature of the subsequent hot, radiation-dominated phase. 
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1. Introduction 

Measurements of the microwave background temperature made by the WMAP satellite 
over the last decade [0 are consistent with the idea that structure in the universe origi- 
nated as small fiuctuations during a primordial era of accelerated expansion, or 'inflation'. 
Inflation leaves the universe in a cold vacuum state, devoid of radiation and matter. If 
our own universe experienced an inflationary era in its past, we must conclude that some 
process repopulated its sterile vacuum with the abundant quanta we see today. Much 
less is known about this 'reheating' process than is understood about the observationally 
accessible density fluctuation. 

Properties of inflationary fluctuations largely decouple from the details of whatever 
physics drives the accelerating epoch, leaving only weak traces.^ This freedom has en- 
couraged many proposals for the underlying microphysics, nearly all of which are equally 
unconstrained by observation. In an especially interesting proposal, due to Silverstein & 
Tong P] and later elaborated by Alishahiha, Silverstein & Tong 0], inflation occurs while 



""^These traces are presently the subject of intense theoretical and observational effort; see, e.g. Ref. 



a D3-brane moves in a warped background spacetime.^ As we shall discuss in §y, the 
motion of the brane is controlled by an action of Dirac-Born-Infeld ( "DBI" ) type and the 
resulting model is known as DBI inflation. In this paper we study the reheating era which 
must follow such an inflationary phase. Reheating in related brane-world models has been 
studied in Refs. |^ ^, ^, 11, 12|. In contrast to the primordial perturbations, reheating 



depends strongly on many model-specific details |13, 14|. It is a subtle function of nonlinear 



physics, probing information complementary to the linear physics of the density fluctua- 
tion. Nonlinear properties vary widely between competing inflationary models, and there 
are important differences between DBI models and those with canonical kinetic terms. 

Many authors have contributed to the theory of particle production, which is now 
well-developed p^ , ^, 17, 18|.^ In perturbative reheating, the decay rate of inflaton parti- 



cles into other species of matter is calculated directly from an S-matrix. These calculations 
assume any relics of previous decays to be sufficiently diffuse that particle production takes 
place into an effectively pristine vacuum. It was later understood that a Bose enhancement 
of the decay rate may occur for integer-spin species if the products of previous decays accu- 
mulate, leading to a resonant phase of out-of-equilibrium production known as preheating. 



A different effect, so-called tachyonic preheating |^, 24, 25|, is associated with the end of 
inflation in certain models. Unlike conventional preheating this does not merely convert 
one species of particle into another, but (like any tachyonic instability) reflects a preference 
to convert one vacuum state to another by the rapid accumulation and condensation of 
particles. 

These phenomena may occur in any effective field theory of the post-inflationary uni- 
verse, whatever its microphysical origin. The general theory has subsequently been applied 
to concrete models, including some versions of brane inflation, in which it is possible to 



give a more refined interpretation |26, ^. In these models, inflation ends when the moving 
D3-brane becomes close to an antibrane (a D3-brane), allowing open strings stretching be- 
tween the D3-D3 pair to become excited. The open string states include a tachyon which 
induces fragmentation of the D3 brane into DO branes, described on the D3 worldvolume by 
a phase of tachyonic preheating. The DO fragments decay into closed string states, which 
are enumerated by the Kaluza-Klein ("KK") modes of supergravity fields in the warped 
background. These KK modes subsequently decay into the species of matter and radiation 
which populate our universe [p7|l . 

In this paper we study a version of the DBI model in which the end of inflation is 
less dramatic. We work in a model where the D3 comes to rest near an extremity of its 
warped background, referred to as the "tip of the throat." As it settles on its resting-place 
it executes coherent oscillations, inducing growing fluctuations which fold and wrinkle 
its surface. This scenario was suggested by Silverstein &: Tong |^]. Our interest is in 
the qualitative differences which may emerge in concrete models. Should we expect the 



^D-branes are extended objects in spacetinie. For a review of their properties see Johnson [H]; inflation 
in such models has been reviewed by McAllister & Silverstein H. Many detailed properties were discussed 
by Chen §. 

^The theory of preheating has been reviewed by Allahverdi et al. ||19|, Kofman et al. fcoj] and Bassett 
et al. Ell. See also Ref. H. 



reheating temperature to increase or decrease? Is parametric resonance more efficient, 
perhaps due to the extra non-hnearities of the DBI action, or do these same non-hnearities 
conspire to shut off the resonant channels? For this reason we focus on differences with 
the case of canonical kinetic terms and largely ignore similarities. For example, the use 
of perturbation theory in either case will break down at the onset of back-reaction, when 
the population of matter species grows to the point where it cannot be ignored. This 
affects the canonical and DBI cases equally, and we restrict our comparison to times before 
back-reaction becomes significant. 

This paper is organized as follows. In §|2| we review the dynamics of infiation in DBI 
models and discuss a toy model in which inflation ends as the brane settles to its final 
position in the warped background. This model is comparable to chaotic inflation in a 
theory with canonical kinetic terms. In §^ we discuss parametric resonance, making a 
comparison with canonical chaotic inflation. In §^ we study the stage of reheating and 
estimate the reheating temperature. We conclude in §|5[ 

Throughout this paper, we work in units where h = c = 1 and set the reduced Planck 
mass to unity, Mp = (SttG)^^'^ = 1. The metric signature is taken to be (— , +, +, +). 

2. Dirac— Born— Infeld inflation 



In §2.1 we review inflation using a Dirac-Born-Infeld Lagrangian and study the oscillating 



era which occurs while the brane is settling near the tip of throat. In §2.2 we discuss 
possible couplings of the inflaton to matter. 

2.1 Review of DBI Cosmology 

In the scenario discussed by Silverstein & Tong [^ , one or more D3-branes move in a warped 
throat, taken to be an elongated region of a six-dimensional compact manifold with metric 
dsg. Transverse distance within the throat is measured by a parameter r, and the entire 
ten-dimensional geometry can be written schematically 

ds?o = h-^/\T)rj^,dx^'dx'' + h^^\T)dsl (2.1) 

where /i(r) is a warp factor. As r — )• oo there is an ultraviolet regime where the throat is 
glued onto to the rest of the dsg. In this region the throat corresponds to a cone over an 
Einstein manifold. In the infrared regime where r — )■ 0, corresponding to the "tip" of the 
throat, fluxes wrapped along cycles of the conifold smooth the conical singularity with an 
S^ cap. 

A low energy description of a D3 brane travelling in this warped throat is given by 

S = ^Jd^x./^[R + 2P{(l),X)Y (2.2) 

where Qf^^ is the pull-back of Eq. ( ^.1| ) to the brane worldvolume, and R is the Ricci scalar 
constructed from g^i,. The worldvolume scalar <j) is a collective coordinate which measures 



the radial position of the brane within the throat, and X = —g^'^dfj_(j)du4>- The action 
describing evolution of (j) satisfies 

P{4>,X) = -fi{(t>) [Vl - fM)X - l] - V{(t)), (2.3) 

in which /i and /2 are determined by the warped geometry and y is a potential which 
describes how the brane is attracted towards the bottom of the throat. By choosing coordi- 
nates we can arrange that /i = /2~ = /^^- In these coordinates f~^ measures the tension 
of the D3-brane, redshifted by the warp factor, giving / canonical dimension [GeV~ ]. The 
proper tension T3 is set in terms of the inverse string tension 27ra' and the string coupling, 
Qs, by Tg-i = {2T,f{a'fgs. After redshifting, / = h/T^. 

If the throat is exactly anti de Sitter, the result is a static BPS state in which the 
Coulomb potential is quartic. Conformal invariance forbids the existence of a mass term. 
If the throat is cut off then conformal invariance is broken, and we expect a mass term 
to be generated which receives contributions from moduli stabilization and bulk fiuxes 



pq , 28, 29, ^, 31, |3^, ^. In a general compactification, a spectrum of terms will be 
generated including a (j)^'"^ term which may lead to a phase of 'ledge' or inflexion-point 
inflation. We assume this term is absent, or forbidden by a symmetry. Therefore, for 
small departures from the equilibrium point, the potential attracting the brane towards 
the infrared tip of the throat is quadratic Q , 

V{<^) = \rr?<l?. (2.4) 

D-brane inflation with this potential is sometimes referred to as the 'UV model'. A related 
small-fleld (TR') model due to Chen [^ takes V = Vq — m^cj)^ /2 with Vq dominant, but in 
this paper we do not consider this possibility. In DBI models there is less need to tune the 
contributions to m^ to avoid an rj problem [p^], since inflation can be supported on steeper 
potentials than would be allowed with canonical kinetic terms. 

As the brane drifts towards the infrared tip of the throat, slowly decreases to zero. 
In this limit (jP' -^1 and the square root can be expanded in powers of X, reproducing the 
action for a canonically normalized scalar field with potential V . On the other hand, if 
the brane is moving relativistically then X ~ 1 and the square root in Eq. ( |2.3| ) cannot be 
reduced.'^ To measure the relative importance of these nonlinear terms it is conventional 
to define a Lorentz factor for the brane, 



l-'^^Jl-f{<P)<t^^. (2.5) 

Unitarity of Eq. ( p.2| ) requires 7 to be real, and therefore ^ fcj)^ < 1 [^8| , |3|, which gives 
a limiting velocity for (j). For a slowly drifting brane we have 7 ~ 1, whereas relativistic 



*In this limit we retain an infinite set of non-renormalizable terms involving powers of {d4>Y with 
coefficients in pre-determined ratios, but discard all operators of the form d^<j) for n > 1. This is nevertheless 
sensible because yjl — fX is the unique operator containing only dtf) which is compatible with higher- 
dimensional isometries of the DBI model, and therefore satisfies a form of the non-renormalization theorem 



motion corresponds to 7 ^ 1. The equations of motion for the inflaton are 



+ ^H^-i2 + oT^' + is(l2+y']=^^ (2-6) 



where a prime ' denotes a derivative with respect to (j). In the non-relativistic hmit / — s- 



and 7—7-1, and Eq. (2.6) reduces to the Klein-Gordon equation for a canonically normal- 
ized scalar field. 

Maps of the microwave background temperature anisotropy place stringent constraints 
on 7, which can source a large signal in the equilateral mode of the bispectrum of order 
/nl 0.3272 ||. An analysis of WMAP 5-year data yields /nl = 155 ± 140 in this 



channel |3£], where the quoted error is la; at 95% confidence the requirement is roughly 
7 < 18. In this paper we will assume this upper limit will continue to be eroded, perhaps 
ultimately constraining 7 to be rather small, so that it is natural to imagine the model is 
close to its non-relativitic limit. Accordingly we suppose that inflation is approximately of 
the slow-roll variety. In this limit the evolution of (j) is close to a canonical model, up to 
small corrections controlled by the deviation of 7 from unity. We will solve Eqs. ( |2.6| )~(^^) 
perturbatively in powers of these small deviations. 

Inflation ends near the tip of the throat, where the S^ cap smoothly cuts off further 



decrease of the warp factor pO| , 41 1. In this region there is no redshifting associated with 
motion down the throat, and the brane tension /^^ becomes approximately constant. 
Silverstein & Tong [^ suggested that the relevant physics of the capped throat could be 
captured by approximating the tension as 

/('A) = ... ^ ..2 > (2-7) 



A2 



+ /^2) 



where A is the 't Hooft coupling and /i is the energy scale of the infrared cut-off. We 
have chosen the origin of (j) so that the final resting place of the brane can be taken to 
be (/) = 0. The warp factor corresponding to Eq. ( |2.7D is not itself a solution of Einstein's 
equations, but is qualitatively reasonable. When (j) <^ fi, this becomes independent of the 
brane position (j), 

/^^. (2.8) 

We will discuss natural values for A and /i in §§. 

We restrict our attention to small oscillations around this position, contained entirely 
within the cut-off region where /"^ is approximately constant, and suppose that m ^ H. 
Therefore many oscillations take place in one Hubble time and terms in Eq. ( |2.6| ) arising 
from Hubble friction can be ignored. We will justify these statements below. Making all 
these approximations and taking the non-relativistic limit where 7 = 1, coherent oscillations 
of the brane occur with constant amplitude A, 

0o(t) ^^cosmi. (2.9) 



Validity of perturbative approximation. Our analysis depends on the approxima- 
tions described above, and will provide a useful description only when perturbative rela- 
tivistic corrections dominate those from other sources. 

When is this likely to occur? After integrating out fluxes and other UV perturbations, 
we are left with a low energy theory controlled by an effective potential for the D-brane 
p^ , [3^ , of the form V = J2i Cirn^~^(f) \ where the Cj are expected to be 0(1) in a generic 
compactification. The mass m is an ultraviolet scale which determines where the throat is 
glued onto the main bulk of the compactification manifold. The leading term corresponds 
to A = 3/2. Assuming some symmetry forces this to be absent, the next A = 2 term will 
typically dominate higher contributions when (f) <^ m. 

The infrared cutoff scale, fi, must satisfy fj, < m. Therefore, provided we take (j) < fi, 
in order that we are safely within a region where the warp factor is constant, we can 
expect corrections from both features in the throat geometry and higher-order terms in the 
D-brane potential to be negligible. 

Now consider our neglect of Hubble drag in Eq. (^Is|). In practice, Hubble drag would 
cause A to decrease slowly, roughly like (mt)^^ [42|. Suppose Hf and Hi are two values of 



the Hubble parameter separated by A^ oscillations. Since the brane worldvolume experi- 
ences a matter-dominated era while the coherent oscillations are ongoing, these values are 
related by 



-±r^l + AN—, (2.10) 

Hf m 

where ^ is a constant of order unity. The corresponding oscillation amplitudes are 



in which A' is a different constant, also of order unity. Unless A^ ^ 1, our assumption 
that m^ H is sufficient to ensure Hi ~ Hf. In addition, Eq. ( |2.11| ) shows that Ai <^ Mp 
is sufficient to ensure Ai ~ Af. In view of the requirement A <^ m required to ensure 
the quadratic potential is an adequate approximation, where m is typically of order the 
infrared cutoff scale, this is amply satisfied. Therefore, in what follows, we will restrict our 
attention to A^ < 10^ oscillations and neglect the time dependence of H and A. 

Relativistic corrections to background solution. Corrections to the non-relativistic 
motion described by ( p.9D will be perturbative provided |/(/'ol ^ 1. Introducing a dimen- 
sionless positive quantity e, satisfying 

e = fm'^A^, (2.12) 

this perturbative condition can be rewritten e ^ 1. We must nevertheless work in a regime 
where e > (j)/m or /u/m, in order that the small relativistic corrections we obtain are not 
dominated by others which we do not calculate. We should trust our analysis for perhaps 
e ~ 10~^ to 10~^, but not for significantly larger or smaller e. 

To go further, we must determine the leading 0(e) corrections to (^.9|). Note that e 
is not the conventional slow-roll parameter, given by —H/H^, and these corrections are 



not related to the slow-roll approximation. The governing equation can be obtained by 
dropping Hubble friction terms in Eq. (|2.6|), yielding 



+ 



-I 



dV 



(2.13) 



It is tempting to imagine that Eq. ( |2.13 ) can be solved perturbatively, writing (/) = (/)q + 54> 
with 6(j) of order 0(e). However, Eq. (|2.13 ) is rather nonlinear and also induces corrections 
to the fundamental frequency m. To capture these effects, consider a WKB-type solution 
of the form </> = Acos9{t) and solve for the phase 9{t) to order 0(e). To match smoothly 
with ( ^.91) in the limit e — )- 0, we must have 9{t) = mt + 60{t). From Eq. ( 2.13 ) we conclude 
that 69 (t) must satisfy 



i?77 r 

69 + 2m(cot mt)69 H sin mt cos mt = 0. 

Dropping an unwanted divergent contribution,^ the solution is 



(b(t) = A cos <i mt (1 I H sin 2mt 

^^ ^ ^ * 16/32 



(2.14) 



(2.15) 



Corrections to the amplitude are also calculable in principle, but need not be computed if 
we are ignoring the time dependence of A. 

In comparison with the canonical solution, the phase shows a secular drift proportional 
to emt. Although the calculation remains under perturbative control, this drift generates 
order unity deviations from the non-relativistic background after 1/e oscillations. There is 
also a superposed phase oscillation, proportional to |e sin 2mt\ <C 1, which always generates 



small corrections. To make contact with the canonical solution, (2.15) can be written in 
the alternate form 



cPit) 



AU-- sin2 
16 



, -, 3e 

mt 1 

16 



, . 3e 

cos mt 1 

16 



(2.16) 



which is equivalent to ( 2.15 ) up to 0(e). Eq. ( 2.16| ) exhibits the first perturbative relativistic 
corrections to the motion of the brane: first, a frequency shift in the background motion (J)q, 
corresponding to m i— )• ?n, (1 — 3e/16); and second, a small periodic fluctuation superposed 
over the background oscillations. 



2.2 Coupling to matter fields 

The reheating phase is designed to fill the universe with Standard Model degrees of free- 



dom. Pauli's exclusion principle usually prevents efficient preheating into fermions [43|, 
so we suppose the inflaton <j) couples to bosonic modes which are taken to have spin zero. 
The resulting scenario is frequently adopted in the literature, and gives an acceptable phe- 
nomenology subsequently decay to fermions is efficient. We treat the couplings between 



^This mode is proportional to cot mi, which is divergent at f = 0. It is projected out by the boundary 
conditions, and therefore does not lead to an instability. 



the inflaton and the matter modes phenomenologically, denoting the matter field x ^'^'^ 
taking it to be described by the low energy effective Lagrangian 

C 5 -\d^xd^x - \mlx^ - \g^cf'x\ (2.17) 

where (7 is a dimensionless coupling and m-^ is a bare mass. The symbol '5' is used to 
indicate that the Lagrangian contains these contributions among others. We do not claim 
that such interactions necessarily arise in string models, although we hope the resulting 
phenomenology will be representative of those that do. The equation of motion for x is 

X + 'iHx-^V\+{ml+g^ <!)-') X = Q. (2.18) 

where V^ is the spatial coordinate Laplacian. Matter production in this model is studied 



in S3. 2 



In §3.3 we will consider an alternative model, in which the matter field x is coupled to 



a symmetry-breaking potential of the form 

1/(0, x) = \rr? (0 - cjf + ^<7'0'x'- (2.19) 

This causes the brane to be attracted to a radial position for which (p = a, representing a 
toy model of compactifications in which the mobile D-branes are brought to rest at a finite 
point in the throat. 

3. Parametric resonance and particle production 

As the brane settles towards its final resting place, its coherent oscillations lead to copious 
non-perturbative particle production. We briefly review this process, with the aim of fixing 
our notation. 

To quantize any field, such as the inflaton (p or the matter field Xi we pick two inde- 
pendent solutions to the mode equation, 

Xk + 3Hxk + {k'' + ml+g^^^)xk = 0, (3.1) 

where k is the comoving wavenumber, and to be concrete we have chosen to frame our 
discussion for x although an equivalent construction can be given for cj). Label these 
solutions x'k ^^"^ Xfc*- The Heisenberg field corresponding to x can be written 

X{t, x) = / (03 {xfal + xf a_,) e'"^-, (3.2) 

in which we have assumed that the normalisation of {Xk^Xk*} ^^^ been adjusted so that 
Ofc and al, obey the usual creation-annihilation algebra, 

[afc,4,] = (27r)3<5(k-k'). (3.3) 

What is the occupation number n^ for momenta of magnitude k? The energy density, 
Ek, associated with x-Quanta of wavenumber k satisfies 2Ek = IXfc P + '^llXfc Pj where 



o;^ = fe^ + g^cfP'. Accounting for subtraction of the energy density associated with vacuum 
fluctuations, when the occupation number is large it must be given to a good approximation 

by 

n. « ^ = ^ f ^ + \xt?] - \. (3.4) 



ujk 2 



bJ 



It follows that the occupation number can be determined from knowledge of the mode 
functions x'k ^'^d their derivatives; when \x'k\ o'^ \Xk\ ^^^ large, we expect large occupation 
numbers. In what follows we shall not need to make use of the quantum field x? which 
allows us to drop the label 'cl' on each mode function. In particular, if Xk experiences 
exponential growth — implying that the production of x particles is stimulated by the pres- 
ence of existing x-quanta — then Eq. (3.4) shows we can expect n^ to attain macroscopic 



values very rapidly. This is the phenomenon of preheating, discussed in §||. 

To determine when Xk grows, it is often useful to translate the mode equation for Xk 
into a version of Mathieu's equation or its generalizations. The classical Mathieu equation 



is m, |5| 



J 2 

— ^ + {Ak - 2q cos 2z)xk = 0, (3.5) 

dz^ 

where we have introduced a dimensionless time z to be defined below. In the following 

sections we give explicit expressions for A^ and q applicable to both inflaton and matter 

fluctuations. For certain combinations of these parameters the solutions of Mathieu's 

equation exhibit exponential growth. If q > A^, the effective mass can become temporarily 

zero. In this case it is energetically inexpensive to copiously produce x-quanta, referred to 

as broad resonance. Alternatively, if q <^ A^ then the effective mass never passes through 

zero, and the combinations which yield exponential growth occur only in narrow bands of 

momentum space. This scenario is known as narrow resonance. These bands are centred 

on momenta for which A^ = t"^, where £ G N, and roughly correspond to values of k for 

which \Ak — f'\ < q^. If 1^1 <C 1 then at large £ these bands become increasingly narrow. 

3.1 Inflaton production 

The foregoing analysis applies equally to inflaton and matter fluctuations. First, consider 
the inflaton mode of wavenumber k. We set (j) equal to its background value, denoted 
(/?, which is determined up to first-order relativistic corrections by ( ^.161) . At the onset 
of oscillations it is reasonable to assume that x is in its vacuum state x = 0. Until the 
accumulation of x-Quanta causes the background x field to grow, the infiaton fiuctuation 
6(j)k evolves according to the Klein-Gordon equation, to linear order in e. 



3:^^^54+f4 + ^^)<^0. = O. (3.6) 



7 V7' 



We expand uniformly to leading order in fip^ and make the change of variables 54>k = CUk^ 
choosing ^ to eliminate the friction term in ( |3.6D , 



ding 3 2- /o 7^ 

-dT = 2^"^ ^^- (^-^^ 



The solution to this differential equation is 

e = eoexp l^/m^ [^\t) - ^\to)] I , (3.8) 

where (, = ^q at t = tQ. This is an effective friction, or 'braking,' which could be interpreted 
as a perturbative manifestation of the D-cceleration mechanism [^ . It is a novel effect in 
the DBI scenario, with no analogue for canonical kinetic terms. Accordingly, the braking 
disappears in the limit e — )• 0. The remaining time dependence is carried by y^, which 
satisfies a Mathieu equation with Af^. and q determined, to leading order in e, by 

/c2 \ / 3e\ 3e , 3e , , 



The associated dimensionless time z, which was used to write Eq. (|3.5| ), satisfies z = 
mt{l — 3e/16). Evidently q <^ A/., giving resonance only in narrow bands. The i = 1 band 
is centred at momentum ki, given by 

m 3e 

According to Floquet's theory each mode behaves like e^''^. In the i = 1 band, the Floquet 
exponent //^ has positive real solutions which correspond to exponential growth [^] 

Ml = (|)'-(^f-l)'. (3.11) 

This exponent is maximal for the central value k = ki and decreases for larger or smaller 
k, forming a band of width Ai, 

At the centre of the band, the Floquet exponent takes the value fii 

m = max Hk = — . (3.13) 

e=l band lb 

Up to this point, we have written our formulae in terms of the mass, m, and coupling, 
g, which appear in the Lagrangian. However, these quantities only parametrize the theory. 
To make a meaningful comparison between the DBI model and the case of canonical kinetic 
terms, we must express our answer in terms of a measurable mass and coupling. We denote 
these measurable quantities m* and 5*, to distinguish them from the Lagrangian param- 
eters. In the DBI theory the Lagrangian parameters have no direct physical significance, 
whereas in the case of canonical kinetic terms they coincide with the measurable mass and 
coupling at tree-level. 

We determine the relationships m^(m,) and gi,{g) in Appendix |A[ Collecting the results 
of Eqs. ( ^.8D and ( p. 11 ), and expressing the answer in terms of m,^ and g^,, we find that up 



10 



to an overall normalization near the centre of the resonant band, the inflaton modes grow 
like 



d(f>k ~ exp < — cos 



+ ^m*t } . (3.14) 



7e 
2m,t I 1 - - 

The undetermined normalization absorbs the constant ^q iii ( |3.8| ). 

At early times, the DBI friction and the resonant growth compete, and the argument of 
the exponential in Eq. ( |3.14| ) can even become negative. After a few oscillations, however, 
the cosine which represents the friction effect becomes negligible in comparison with the 
resonant term. 

3.2 Matter particle production 

Our ultimate interest lies with the production of matter particles. Consider the poten- 



tial ( 2.17 ). For any particle species which gains a mass at energies much less than the 
inflationary scale it will typically be a good approximation to assume rriy. is negligible in 
comparison with (70 unless the amplitude of oscillations is extremely small. Therefore each 
X-mode receives an effective mass of order k"^ + g'^(jp'. As in §^?^ we assume that preheating 
begins in the vacuum x = 0) which will be valid until the accumulating ^-quanta backreact 
on the zero mode. In this vacuum, small fluctuations grow according to a Hill equation of 
type 

^Xk + (^0 + 26*2 cos 2z + 26^4 cos Az) 5xk = 0, (3.15) 

where the parameters ^O; ^2 and 64 obey 

e / 3e\ g^A^ / 9e 

9'^'^l + ^] (3.16) 



4m^ 

12 

eg', 



_3_^ 

128 m"^ 



and (as in §3.1) a prime ' denotes a derivative with respect to ^ = mt{l — 3e/16). The 
solutions of Hill's equation grow and decay with characteristic Floquet exponent fi^ in a 
way entirely analogous to the solutions of Mathieu's equation. The exponentially growing 
solutions organize themselves into bands, characterized by integers i. Using the arguments 
of Floquet theory one can obtain an estimate of Hk for the i = 1 band (see Appendix y) , 

We believe this expression has not previously appeared in the preheating literature. In the 
region where our approximations are valid, we have A <^ m. Therefore, the terms involving 
9l and ^f are negligible. Under these conditions the Floquet exponent can be simplified, 
yielding 

^i-(|)'-(^r-i)' • (3.18) 
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The centre of the 1 = 1 band hes at k = ki, for which 9o{ki) = 1, 



2 /12 



k 



m 



g'A 



3e 

-r I fn 



2 /12 



g'A 



(3.19) 



Collecting terms as before, it follows that the unstable matter modes grow at a rate deter- 
mined by 



5xk ~ exp 



^^k\ 



m^,t 1 



3e 

16 



exp 



m*t 1 + — 

m* / V 16 



(3.20) 



This displays a slight enhancement compared to the equivalent canonical model. 
3.3 Preheating with a symmetry breaking term 



Now consider the symmetry-breaking potential given in Eq. ( 2.19 ). It is convenient to 
redefine (p so that it describes excitations around the true vacuum at = cr. After making 
the shift <J3 ^ <j) + a, the potential can be written 



vicj), x) = ^rn^cP^ + ^g^^P^^ + g'^^'Px^ + ^a^'^^x^- 



(3.21) 



Excitations of x around this minimum acquire an effective mass rn^^ {(p) = m?' + g^X^ ■ 

Making the same assumption that preheating begins while x ~ 0, the inflaton modes 
will behave as in Eq. ( 3.14| ) until the growth of x occupation numbers causes the x ^^ro 
mode to grow. Until that time, the x modes obey 



^Xk + {k^ + 5^' + 25V(/. + 5^2) 8xk = 0. 



(3.22) 



We will suppose the brane executes oscillations around (j) = a with amplitude ^ <C o", so 
that it can reasonably be said to have settled at its final position. Assuming that a <^ rn/g, 
and again neglecting any bare mass for %, fluctuations in the matter fleld will be controlled 
by a Hill equation, 

^Xk + [^0 + 26*2 cos 2z + 206 cos 6z] 5xk - 0, (3.23) 

in which 2z = mt{l — 3e/16) and the parameters Oq, 02 and 84 obey 



.„ = i,(p+»v)(i + f) 



gagAf 21e 
m m \ 64 
3e ga gA 

h — • 

16 771 m 



(3.24) 



To determine the growth of each x-mode we require an estimate of the Floquet exponent 
associated with Eq. ( |3.23| ) , which is again derived in Appendix |^. For the i = 1 band, it 
is given by 



fJ-k ^ 



92 , 



- 
32 ^ I 2 



,1/2 
'0 



(3.25) 
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As above, the terms involving 9^ and ^f can be neglected. The result is equivalent to a 
Mathieu equation with coefficients A^ and q satisfying 



Ak 



m^ 



{k^+g^a^)(l + j] and q 



92 



.gfy 9^ 



m m 



1 + 



21e 



<^Ak 



(3.26) 



Since q is much smaller than A^, we expect resonance to be of the narrow type if it occurs 
at all. The first instability band occurs when A^ = 1 and is centred at wavenumber k = ki, 
where 

which is shifted by 0(e) in comparison with the canonical case. Eqs. ( [3.19 ) and ( |3.27| ) show 
that the i = 1 resonance injects energy into a narrow band of modes with momenta k ^ m. 
Kofman, Linde & Starobinksy interpreted the narrowness of this resonance as multiple, 
rapid copies of a decay process in which a single inflaton particle of mass ~ m decays into 
two x-particles with opposite momenta of magnitude ~ ki |l^. The width of the resonant 
band satisfies 

Afc,:.2m^^fl + ^V (3.28) 

m m \ b4: / 



Using the appropriate limit of Eq. ( 3.25 ) to estimate /ifc, or applying ( 3.11 ), we conclude 
that there is a small enhancement of preheating efficiency owing to the relativistic DBI 
correction. Near the centre of the i = 1 band, the ^-modes grow like 



^Xk 



exp(,/ifc| 



. z) ~ exp 



m^i 1 H 



(3.29) 



in which the growth rate is augmented with respect to the non-relativistic or canonical 
cases by an 0(e) term. 

3.4 Efficiency of narrow-resonance preheating 

The resonances we have discussed above are narrow, and therefore occupation numbers 
grow only for a limited set of momenta. This process is much less efficient than broad 
resonance, where a large range of wavenumbers experience growing occupation numbers. 
Nevertheless, because growth is exponential within the unstable bands, preheating by nar- 
row resonance may still lead to an acceptably rapid conversion of inflaton modes into 
matter species. Gravitational redshift will cause modes of fixed comoving wavenumber to 
drift slowly through each unstable band from ultraviolet to infrared, and significant particle 
production will occur only when modes remain within an unstable band for sufficiently long 
that multiple e- foldings of stimulated emission occur [18|. If wavenumbers move through 
the unstable bands too rapidly then little stimulated emission occurs, and the resonance is 
washed out. 

The arguments of Floquet's theory show that in each resonant band, the inflaton and 
matter fluctuations of wavenumber k grow like e^*^. Let us focus on the broadest instability 
band, which corresponds to i = 1 and will generally dominate the particle production. This 
band has width ~ q, so each wavenumber remains within the unstable band for a time of 
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order q/H. The resonance will be stable against Hubble expansion if [^i^z/t) q/H > 1. It 
follows that preheating will drain energy from the inflaton provided 



2i? / 3e , , 

For the matter field, preheating will operate in the cases (7 = and o" 7^ whenever 



or 



2H ^ 1 fg^,A\ ( 13e\ 



m^ ~ V m^ y V "^* y V 64 y ' 




apply, respectively. 



4. Perturbative reheating 

While preheating is occurring it remains possible for inflaton particles to decay into other 
species via perturbative processes |Q. Indeed, these processes are necessary to thermalize 
the post-inflationary universe and bring reheating to an end. 

In the theories we have been considering, such decays are allowed by the symmetry- 
breaking potential studied in §p^, where there is an effective interaction g^acjix^ in the 



true vacuum which permits the decay (j) — t- xX- The simple potential ( 2.17 ) allows only 
two-body scattering and annihilation processes in vacuum. However, in the presence of an 
oscillating background field we can take the decay rate to be similar to that produced by 
the g'^(T(j)x interaction with o" ~ ^. If the inflaton has direct couplings to fermionic species 
then perturbative production can proceed via Yukawa interactions such as h(j)ilj'ip, where 
'(/' is a Dirac fermion and Tp = iip^'j^ its adjoint. In this case we must suppose that any 
interactions in which the inflaton participates are not sufficient to spoil the possibility of 
successful inflation. 

In this section we will continue to assume that the inflaton mass is much larger than 
the bare mass of %, or the mass of any fermionic species where such couplings exist. Under 
these conditions, the tree- level decay rates for cj) — )• xX and cj) — )• -ip^ can be calculated 
straightforwardly. The necessary expressions were given by Kofman et al.^ |18] 



r(0^XX) = |^, and T{<P^i,iP) = -^, (4.1) 



where /i^ is the measured Yukawa coupling, determined by matching Eq. (O) to observa- 
tion. Can perturbative decays give a contribution which is competitive with preheating? In 
one oscillation, perturbative decays populate the x field at a rate T/rri-i,. In the same time 
interval preheating increases the population at a rate ^k- Therefore perturbative decays 
dominate if 

r > m.fik. (4.2) 



^Modifications to these rates and cross-sections in a DBI model are studied in Appendix M. 
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Let us suppose that perturbative decays dominate the production of matter species. These 
processes are relatively slow and it is reasonable to assume that the decay products ther- 
malize. After thermalization has gone to completion we can expect the resulting radiation 
dominated universe to have a temperature Tr, defined by 



Tr ~ 0.2yfMp, (4.3) 

where F is the rate of the dominant decay channel. 

Let us briefly recapitulate the role of perturbative processes in canonical inflation. 
Inflaton decays (j) — >• XX can become the dominant channel for inflaton decay only in 
the case of symmetry breaking, for which a 7^ |]l^]. These decays begin to dominate 
when the amplitude of inflaton oscillations has dropped to ^ ~ g'^a/Sir. If we suppose 
that the curvature perturbation synthesized during the inflationary phase, denoted C, was 
dominated by fluctuations in the inflaton, then we must conclude m ~ 10~^Mp in order 
to match observation. This assumption is minimal but not mandatory; if C receives a 
dominant contribution from other sources, a wider range of masses m can be tolerated. 
Supposing further that the coupling g^, has a natural value 5^ ~ 10~^, we find a reheating 
temperature of order the symmetry breaking scale 

Tr ~ a. (4.4) 



Alternatively, it may happen that the dominant decay is into fermions. If so, Eq. ( [4.2| ) 
implies that the Yukawa coupling must satisfy 

hl>mJ^)(^). (4.5) 

Suitable values for the inflaton mass m^ and couplings {g^,hi,} can be found so that this 
condition applies, but these choices are model-dependent. As an example, if we assume 
/i^ ~ 10^^ and m^ ~ lO^^Mp, we find the universe reheats to a temperature 

Tr ~ 1.3 X lO^^A/p ^ 3.12 X 10^^ GeV. (4.6) 

4.1 Reheating temperatures 

The above analysis can be applied at once to the production mechanisms discussed in 
§§p?l-3.3. Let us suppose that cr 7^ 0, so that we are working with the symmetry-breaking 



potential. Assuming a typical coupling g^, ~ 1, Eq. (4.2) implies that — )• xX decay 



dominates resonant production when the amplitude becomes of order 

^ < 4 X 10-V (1 - ^) (4.7) 

or smaller. The resulting reheating temperature receives relativistic corrections in compar- 
ison with the temperature achieved using canonical kinetic terms: 



Tr ~ 0.13CT 1 - - . (4. 
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On the other hand, decays into fermions will become predominant {(p — )■ ipip) if the associ- 
ated Yukawa coupling is sufficiently large 

''•^-(S5)(^)0-f)- 

This is compatible with the conditions gA/m <^ ga/m ^ 1. In this case, we find that, for 
a representative value of the coupling, say /i^ ~ 0.1 

Tr ~ 3.12 X 10^^ ( 1 - ^ ) GeV. (4.10) 



5. Conclusions 

We have studied mechanisms by which the universe can be repopulated with matter species 
following an era of D-brane inflation. We work in the limit where the motion of the brane 
is at most perturbatively relativistic, which may ultimately be required by observation if 
/nl measured on equilateral conflgurations of the bispectrum is nonnegative. 

Our analysis is valid in a regime where the amplitude of oscillations is significantly 
smaller than whatever infrared scale caps the throat in which the D-brane moves. Indeed, 
because we take the D-brane to be at most perturbatively relativistic, inflation is of the 
slow-roll type. We assume inflation ends as the D-brane gently coasts to its resting position, 
dissipating its remaining energy by executing small amplitude oscillations. The amplitude 
and phase of these oscillations are mildly perturbated in comparison with the oscillations 
of canonically normalized scalar reheating the universe after a phase of four-dimensional 
inflation. 

We encounter some novel effects which are not present in a theory with canonical kinetic 
terms. In canonical inflation there is no resonant production of inflaton particles during the 
reheating era. By contrast, in the perturbatively-relativistic limit of DBI inflation, non- 
linearities arising from the Dirac-Born-Infeld action cause a weak resonance which reaches 
an 0(1) effect only after ~ 1/e oscillations. Excessive production of inflaton particles 
is typically problematic, but since preheating typically concludes after < 10 oscillations 
this weak resonance is unlikely to be fatal for the DBI model. In the ultra-relativistic 
limit where 7 3> 1, the potential becomes unimportant and presumably resonance will be 
strongly suppressed. 

Another novel effect is a 'DBI friction,' which appears in the equation of motion for 
the inflaton. This slightly damps the preheating effect. Beyond perturbation theory, it is 
possible this effect represents the decreasing significance of the potential at large 7. 

For a given set of measured masses and couplings which characterize the properties of 
the inflaton and the matter fleld into which it decays, we flnd that final reheating tempera- 
tures (assuming they are achieved by perturbative decays) are typically smaller than in the 
equivalent theory with canonical kinetic terms. This may be beneficial in concrete mod- 
els, where lower reheating temperatures allow problems associated with overclosure of the 
universe by moduli or gravitinos to be ameliorated. The most stringent constraint on DBI 
models comes from observations of the microwave background bispectrum, which presently 
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require 7 < 18. If our results are representative, and the reheating temperature falls as 7 
increases, a constraint on 7 may also emerge from demanding that the nucleosynthesis era 
is uninterrupted, which requires Trh ^ MeV. 
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A. Canonically normalized theory 

Consider a general higher-derivative Lagrangian of the form P{(j),X), given by Eq. (p.2|). 
We allow small excitations 6(p around the background field (pQ. Assuming derivatives of (po 
are negligible there is a corresponding perturbation in the kinetic energy, X = 5(fP' — {d6(j))'^. 
To quadratic order in S4> but including all orders in other fields, the field propagates in a 
locally Minkowski region according to the action 



S= I d^x iPo + ^54) + 



dP ^ . dP 
dX 



6<p'-{d5cl>r\+^^5cl>' + --- \, (A.l) 



where Pq = P[X = 0, (po) and '• • • ' denotes terms of higher order in 5(j) and its derivatives 



which have been omitted. Eq. (AJ) will give a reliable description of scattering and 
decay processes involving inflaton particles provided derivatives of the background field 
are negligible over the relevant time and distance scales. In this limit we can neglect the 
curvature of spacetime and work in a locally flat region, which should be chosen to be 
somewhat larger than the interaction region. Eq. ( A.lj) applies within this patch. S-matrix 



elements obtained from it can be expected to describe local scattering probabilities. On the 
other hand, if the background field (p varies significantly over the timescale of the scattering 
event then the definition of asymptotic in- and oitt-states becomes more complicated, and 
the concept of well-defined particles with an associated S-matrix may be invalidated.^ 

Canonical normalization. If the vacuum is locally stable on timescales comparable 
with decay processes, then (dP/dcp) must be negligible when evaluated there. We must nev- 
ertheless retain perturbative excitations contained in this term, which include the fermion 
Yukawa couplings. After canonical normalization of 6(p, the interaction terms in the La- 
grangian arise from 

where a comma denotes a partial derivative evaluated in the background. In a cut-off 
throat, the models we have been considering correspond to P = — /~^(1 — fX)^''^ -\- f~^ — 



^For example, this may happen in models where (f>o evolves at a constant rate, such as the ghost inflation 
proposal of Arkani-Hamed et al. wA W8]. However, this model is believed to suffer from certain pathologies. 
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V{cl)) with / constant and V{cl)) chosen to be Eq. ( 2.17 ) or ( 2.19| ). Consider the symmetry- 



breaking potential ( 2.19| ). We only wish to canonicahy normahze 6(j), which represents 



fluctuations around the mean field visible to local observers as particles. Therefore, one 
should shift to canonical normalization only after expanding in powers of fiuctuations 
around the background field. Applying this procedure implies that we should make the 
transformations 

m -^ m, = ——— and g — ^ 0^ = — , (A. 3) 

2Px \APx 

Measurable masses and couplings. In this paper, we study the dynamics of particle 
production in different theories. If we are to make a meaningful comparison, these must 
be theories of the same particles. Therefore our results should be written in terms of 
the masses and couplings measured, for example, by a local observer who can record the 
outcome of scattering events. 

In the case of non-canonical kinetic terms, it is possible that evolution of the back- 
ground field causes the masses and couplings to change in time. In this case, it is the 
masses and couplings measured by a local observer at the time of reheating which are the 
relevant quantities. 

After canonically normalizing the fiuctuation Scj), and working with tree-level scattering 
amplitudes, there are no further rescalings in the LSZ formula. Therefore, at tree-level, the 
coupling (7^ of the canonically normalized field is the coupling which would be measured 
by experiment. Likewise, the mass m^ of the canonically normalized field determines the 
position of the pole in its propagator. It is therefore the quantities m^ and g^, which should 
be compared with the mass and coupling in a theory with canonical kinetic terms. 

If a Yukawa coupling to fermions is present than we should similarly rescale its coupling 
constant, 

B. Floquet exponents from Hill's equation 

In this appendix we derive the Floquet exponent associated with Hill's equation. 

We start with the case when ^2,^4 7^ 0, and ^2™ = 0, for n > 2 and by noting that 
in the edges of the resonance bands Hk = 0; it turns out that for the most important 
instability band and for 02 , 04 < 1 

00 - 1 ± ^2 

which, for 6*4/^2 '^ 1, essentially reduces to 



i^M 

602 



00^1 ±02 . (B.l) 

1/2 

From here, it also follows that, at lowest order, 0^ ~ lib02/2. Furthermore, the instability 



parameter is given by [44] 



sin2(i^,|)=A(O)sin2(0y2|) , (B.2) 
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where A(0) is the discriminant evaluated for fj, = (following from Hill's method); A(0) 
is exclusively determined by the parameters ^O) ^2 and ^4 and is approximately given by 



A(0) :. 1 + ^U. cot ^-^/^ " 



4.9, 



1/2 




'0 



h' ^ e,' 



it can be shown that this, approximately, yields 

2 



A(0) ~ 1 =F 



vr 



-ii)^l(l)^¥ 



(B.3) 



Using ( |B.1| ) and (|B.3| ) , we find that from ( [B.2[) , we may arrive at an approximate expression 
for the instability parameter for this Hill's equation as 



y^fc 



5I/2 
'0 



ii.fm.^ 



(B.4) 



This expression for /x^ is valid for 62,94 <^ 1. 

If we now proceed and consider the case with non- vanishing 9q, ^2 and 9q, we observe 
that the instability parameter may be computed using again using ( |B.2[ ), where now 



vr cot I — 9, 



A(0) ~ 1 + 



^nl/2 



'0 



40, 



1/2 




1 



+ 



9-9, 



(B.5) 



By considering the first instability band, we take ^o — 1 i ^2) so that 
2 ' l-0n ' "^ ' 9-9 ^^' 



g1/2 
^0 



t'2 

1± — 



T^fl±^) and cotf^^^/^ 







-=Ft'92 



(B.6) 



Plugging these into (|B.5|), we find that, up to lowest order 



7r\ 2 
4- 



a(o)-it(t) (T^i + ^ + ^ 



020R ^2 

'2 



Moreover, 



'TT i/2\ , 1 /vr\2 



("J 



/2 



Finally, by substituting into (B.2), we obtain 



fJ-k 



72 \ _ 6'26'g 



T 



32 ^ V 2 



,1/2 



1 



(B.7) 
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